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ABSTRACT 

We performed a timing analysis of the fastest accreting millisecond pulsar IGR 100291+5934 using 
RXTE data taken during the outburst of December 2004. We corrected the arrival times of all the 
events for the orbital (Dopplcr) effects and performed a timing analysis of the resulting phase delays. 
In this way we have the possibility to study, for the first time in this class of sources, the spin-up of 
a millisecond pulsar as a consequence of accretion torques during the X-ray outburst. The accretion 
torque gives us for the first time an independent estimate of the mass accretion rate onto the neutron 
star, which can be compared with the observed X-ray luminosity. We also report a revised value of 
the spin period of the pulsar. 

Subject headings: stars: neutron — stars: magnetic fields — pulsars: general — pulsars: individual: 
IGR J00291-h5934 — X-ray: binaries 



1. INTRODUCTION 

The so-called recycling scenario links two different 
classes of astronomical objects, namely the millisecond 
radio pulsars (usually found in binary systems) and the 
Low Mass X-ray Binaries (hereafter LMXBs), or, at 
least, a subgroup of them. The leading idea of this sce- 
nario is the recycling process itself, during which an old, 
weakly magnetized, slowly spinning neutron star is accel- 
erated by the accretion of matter and angular momentum 
from a (Keplerian) accretion disk down to spin periods 
in the millisecond range. In this way, at the end of the 
accretion phase, the neutron star rotates so fast that it 
is resurrected from the radio pulsar graveyard, allowing 
the radio pulsar phenomenon to occur again despite the 
weakness of the magnetic field. 

Although this scenario was first proposed long ago (see 
e.g. Bhattacharya & van den Heuvel 1991 for a review), 
the most embarassing problem was the absence of co- 
herent pulsations in LMXBs. Only recently, the long 
seeked millisecond coherent oscillations in LMXBs have 
been found, thanks to the capabilities (the right com- 
bination of high temporal resolution and large collect- 
ing area) of the RXTE satellite. In April 1998, a tran- 
sient LMXB, SAX J1808.4-3658, was discovered to har- 
bour a millisecond pulsar (Pspin — 2.5 ms) in a compact 
(Porb — 2 h) binary system (Wijnands & van der Klis 
1998; Chakrabarty & Morgan 1998). We now know seven 
accreting millisecond pulsars (Wijnands 2005; Morgan et 
al. 2005); all of them are X-ray transients in very com- 
pact systems (orbital period between 40 min and 4 h), 
the fastest of which (Pspi„ ^ 1.7 ms), IGR J00291-I-5934, 
has been discovered in December 2004 (Galloway et al. 
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2005, hereafter G05). 

Timing techniques applied to data of various accreting 
millisecond pulsars, spanning the first few days of their 
outbursts, allowed an accurate determination of their 
main orbital parameters. However, only a few attempts 
have been made to determine the spin period derivative 
(Chakrabarty et al. 2003; Galloway et al. 2002). The first 
reported measurement of a spin-up in these sources was 
made for IGR J00291+5934 (Falanga et al. 2005, F05). 

In this paper we apply an accurate timing technique 
to the fastest currently known accreting millisecond pul- 
sar, IGR J0029H-5934, with the aim of constraining the 
predictions of different torque models with good qual- 
ity experimental data. Our results indicate quite clearly 
that a net spin up occurred during the December 2004 
outburst of IGR J00291-I-5934 (see also F05) and that the 
derived torque is in good agreement with that expected 
from matter accreting from a Keplerian disk. 

2. THE TIMING TECHNIQUE 

For a periodic pulsating source, the time of arrival of a 
given pulse at the solar system baricentre is affected by 
three effects which causes temporal delays with respect 
to the predicted ones. These effects are 1) uncertainties 
in the orbital parameters, 2) uncertainty in the spin fre- 
quency and possible secular variations, 3) uncertainties 
in the source position. In standard timing techniques (see 
e.g. Blandford & Teukolsky 1976) the predicted arrival 
time of a given pulse is computed using a first guess of 
the parameters of the system, and the difference between 
the experimental and predicted arrival times, namely 
the residuals, are fitted with a linear multiple regression 
of the differential corrections to the parameters. This 
means that the differential correction to orbital param- 
eters, spin frequency and its derivative, source position 
in the sky, are computed simultaneously. This technique 
has the obvious advantage to give a self-consistent solu- 
tion, where all the correlations in the covariance matrix 
of the system are fully taken into account. However, the 
convergence of the fit is not always guaranteed and - 
especially in the case of long temporal baselines - con- 
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vergcncc to secondary minima might lead to inaccurate 
solutions. 

On the other hand, if the orbital period is much shorter 
than the timescale on which the spin period derivative 
and the source position uncertainty are expected to pro- 
duce a significant effect, it is easy to see that the delays 
in the arrival times produced by the uncertainties in the 
orbital parameters are distinguishable from those caused 
by the other two effects. This is because the former oscil- 
late on the orbital period timescale, while it is expected 
that the latter will follow a secular trend dictated by the 
torques on the accreting neutron star and by the orbital 
motion of Earth. In the following we therefore discuss 
all these effects separately and describe in general the 
timing analysis suitable for these cases. 

The procedure we applied is the following. In order to 
obtain the emission times, tcm, the arrival times of all 
the events, tam were firstly reported to the Solar sys- 
tem barycenter adopting the best estimate of the source 
position in the sky. We note here that the contribution 
of the relative acceleration of the binary system with re- 
spect to the Solar system could in principle produce not 
negligible contribution to the frequency derivative. For a 
source located at the position of IGR J0029i+5934 at a 
distance of few kpc, the most relevant of these effects is 
due to the planar acceleration along the Galactic plane, 
which gives a spin derivative !>Gai ^ 5 x 10~^° (see 
e.g. Damour & Taylor 1991). This is several orders of 
magnitude lower than the spin frequency derivatives we 
found (see below), and therefore in the following these 
kind of effects are not taken into account. Then we cor- 
rected for the delays caused by the binary motion using 
the best estimate of the orbital parameters through the 
(first order approximated) formula: 
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where x = asini/c is the projected scmimajor axis 
in light seconds, m = 27r(tarr — T*)/Porb if' the mean 
anomaly, T* is the time of ascending node passage at the 
beginning of the observation, co is the periastron angle, 
and e the eccentricity. In the following, for simplicity, 
we use t instead of iom- From eq. (1) we compute the 
observed phases as: (j) = v{t — To), where v = i/Pgpin 
and To is the start time of the observation. On the other 
hand, the expected phase variations, 50^,, caused by a 
spin frequency derivative, u, can be computed by a dou- 
ble direct integration: 
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In the simplest case of a constant i), the integration gives 
a parabolic function of time. In general, fitting these ex- 
pected phase variations to the observed ones, we can ob- 
tain an estimate of u{t) and hence important information 
on the torques acting on the accreting neutron star. 

To obtain these information is important to evaluate 
any source of error in the observed phase variations; we 
start discussing the errors induced by the uncertainties 
on the orbital parameters of the binary system. The dif- 
ferential of = ^{t — To), with t given by expression (1), 
with respect to the orbital parameters allows to calcu- 
late the uncertainties in the phases, a^oih, caused by the 
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The uncertainties in the adopted orbital parameters will 

result in a "timing noise" of amplitude cr^orb- These 
should be therefore added in quadrature to the statisti- 
cal uncertainties cr^ stat on the experimentally determined 
phase delays. The resulting uncertainties on the phase 
delays will be: = {al^^^ + (tI,,J^''^.^ 

On the other hand, the uncertainties in the phase de- 
lays, (T^pos, caused by the uncertainties on the estimates 
of the source position in the sky, will produce a sinusoidal 
oscillation on the Earth orbital period. For observation 
times shorter than one year, as it is the case for most 
transient accreting millisecond pulsars, this can cause 
systematic errors on the determination of the neutron 
star period and its derivative, since a series expansion 
of a sinusoid contains a linear and a quadratic term. In 
order to evaluate these effects, let us consider the expres- 
sion of the phase delays induced by the Earth motion for 
a small displacement, 5\ and 5(3, in the position of the 
source in ecliptic coordinates, A and (3 (see e.g. Lyne & 
Graham-Smith 1990): 

A(/)pos = i^oV [ sin(Afo + e) cos 3 5\ — cos(Mo + e) sin fi 6(3 ] 

(4) 

where y = te/c is the distance of the Earth with respect 
to the Solar system barycenter in light seconds, Mq = 

27r(To — Ty)/P0 — A, where is the time of passage 
through the Vernal point, and e = 27r(t — Tq)/Pq « 1 
in our case. 

After some algebraic manipulations, this can be writ- 
ten as: A(/)pos = i^oycr-y sin(Mo -|- e — ^*) w, where a-y 
is the positional error circle, 9* = arctan(tan/3 5(3 /5\), 
and u = [(cos/3 (5A)^ + (sin/3 5(3)^ ^1"^ ja^. Since the true 
source position must lie within the error circle, the follow- 
ing inequalities hold: (cos/3 5\la^Y' ^ 1j {^(^I'^^)'^ ^ !> 
and thus u < {1 + sin^ PY^"^- This means that the uncer- 
tainty in the source position is: 

A(/.pos < I'D 2/fT7 (1 + sin' /?)'/' sin(Mo + €-9*) (5) 
We can expand it in series in the parameter e << 1 in 
order to find the systematic uncertainties induced on the 
linear and quadratic term. Note that, since the values 
of 5X and 5(3 are the differences between the nominal 
and the true (unknown) source position (which can be 
everywhere within the error circle), their ratio is unde- 
termined, and hence 6* can be any value between and 
27r. We have therefore maximized the functions sin0* 
and cos 9* with 1 separately in the linear and quadratic 
terms of the series expansion. The resulting systematic 
error in the linear and the quadratic term of phase de- 
lays evolution versus time, which correspond to the spin 
frequency correction and the spin frequency derivative, 
respectively, are: cTy gyst < '^oy^-y + sin' /3)^/' 2 tt /P® 
and a-£,syst < i^oyo-^ (1 sin' /3)V2 (27r/Pe)'- 

Summarizing, the phase variations caused by /> are ef- 
fectively distinguishable from those induced by the un- 
certainties in the orbital parameters, which result in a 
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"timing noise" of amplitude given by eq. (3). On the 
other hand, the imccrtainty on the source position cannot 
be easily decoupled from the phase variations caused by 
V (particularly for observation times much shorter than 
one year), and therefore results in systematic errors on 
the estimate of the spin frequency and its derivative. 

3. OBSERVATIONS AND DATA ANALYSIS 

IGR J00291+5934 was observed by RXTE between 
2004 December 3 and 21. While the observations be- 
tween December 3 and 6 were already analyzed in G05, 
in this paper we analyze the data between December 7 
and 21 taken from a public ToO. We mainly use data 
from the RXTE Proportional Counter Array (PCA, Ja- 
hoda et al. 1996), which consists of five identical gas-filled 
proportional counter units (PCUs), with a total effective 
area of ^ 6000 cm^, sensitive in the energy range between 
2 and 60 keV. We used data collected in generic Events 
mode, with a time resolution of 125 /iS and 64 energy 
channels. These files were processed and analyzed using 
the FTOOLS v. 5. 3.1. In order to eliminate the Doppler 
effects caiiscd by the Earth and satellite motion, the ar- 
rival times of all the events were converted to barycentric 
dynamical times at the Solar system barycentcr. The 
position adopted for the source was that of the proposed 
radio counterpart (which is compatible with that of the 
proposed optical coimterpart, see Rupcn ct al. 2004; Fox 
& Kulkarni 2004) . For the spectral analysis we also used 
data from the High-Energy X-ray Timing Experiment 
(HEXTE, Rothschild et al. 1998, ~ 20 - 200 keV energy 
range). 

We corrected the arrival times of all the events for the 
delays caused by the binary motion using cq. (1) with 
the orbital parameters given in G05. In order to check 
for the presence of the pulsations during our observation, 
we performed an epoch folding search on each continu- 
ous interval of data (lasting on average 60 min) around 
the spin period given in G05. The pulsation was clearly 
visible up to December 12. After that, owing to poor 
statistics, pulsations could be detected only by folding 

1 day worth of data. No pulsations were detected af- 
ter December 14, in accordance to what reported in G05. 

Using our longer temporal baseline (about 7 days) with 
respect to that in G05 (about 3 days) , we firstly tried to 
increase the accuracy of the orbital period measurement 
using the technique described in Papitto ct al. (2005) and 
successfully applied to SAX J1808. 4-3658. However, no 
significant improvement was found. Adopting the uncer- 
tainties in the estimates of the orbital parameters given 
in G05 in eq. (3), we obtain: cr^orb ^ 0.01, where we 
have maximized sin and cos functions with 1, and used 
t — To < 7 days. Therefore, we expect that the uncertain- 
ties in the orbital parameters will cause a "timing noise" 
not greater than a,po^ x Pgpin ^ 0.02 ms. 

To compute phases of good statistical significance we 
epoch folded each interval of data in which the pulsa- 
tion was significantly detected at the spin period given 
in G05 with respect to the same reference epoch. To, 
corresponding to the beginning of our observations. The 
fractional part of the phase was obtained fitting each 
pulse profile with a sinusoid of fixed period. To c;om- 
pute the associated errors we combined the statistical 
errors derived from the fit, u^stat, with the errors a^orb 

aSCT^ = ((7|,tat+'^|orb)'^^- 



In order to derive the differential correction to the spin 
frequency, IS.vq^ and its derivative, i>o, at the time Tq 
we have to derive a functional form for the time depen- 
dence of the phase delays. We started from the simple 
assumptions briefly summarized below: i) The bolomet- 
ric luminosity L is a good tracer of the mass accretion 
rate M via the relation L = ({GM/R)M, where C < 1, 
and G, M, and R are the gravitational constant and 
the neutron star mass and radius, respectively, ii) The 
matter accretes through a Kcplerian disk truncated at 
the magnctosphcric radius, cx M~", by its inter- 
action with the (dipolar) magnetic field of the neutron 
star. At i?m the matter is forced to corotate with the 
magnetic field of the neutron star and is funneled (at 
least in part) towards the rotating magnetic poles, thus 
causing the pulsed emission. For standard disk accre- 
tion a = 2/7; note that this is very close to the up- 
per limit a = 2.3/7 derived from our data. This upper 
limit on a can be derived noting that the pulsations were 
clearly seen at the beginning of the observations in G05 
at MJD~ 53342.0, and were detected for At ~ 11.5 days 
until MJD~ 53353.5, when the flux reduced by a factor 
of 10. The corresponding expansion in Rj^ must be 
i?MAx/i?MiN = (1 - Ai/^s)-" < Rco/R which, for R = 
10^ cm and m = 1.4, gives a < auAX = 0.328 ~ 2.3/7 
where t|j ~ 12.4 days is the decay time from G05 and 
i?CO is the corotation radius, both defined below. There- 
fore we considered two extreme cases, namely a = 2/7 
and a = 0, since a location of R^ independent of Af 
has been proposed (sec, e.g., Rappaport, Fregeau, and 
Spruit, 2004). iii) The matter accretes onto the neu- 
tron star its specific Keplerian angular momentum at 
i?„i, I = {GMR^y/'^, thus causing a material torque 
Tj^ = ixM. A firm upper limit to this torque is given by 
the condition T^^ < CaxXM, with^ax = (GMi^co)^''^ 
where -Rco = 1-50 x 10* m^/^i^^^/'' is the corotation ra- 
dius (namely the radius at which the Keplerian frequency 
equals u and beyond which accretion is centrifugally in- 
hibited), and m = M/Mq. iv) We do not consider any 
form of threading of the accretion disk by the magnetic 
field of the neutron star (sec e.g. Ghosh & Lamb 1979; 
Wang 1997; Wang 1996; Rappaport, Fregeau, and Spruit 
2004 for a description of the magnetic threading), which 
implies that the only torque acting during accretion is 

'''m ■ 

Under these hypotheses the spin frequency deriva- 
tive is !> = £M/{2wI), where / is the moment of in- 
ertia of the neutron star and we have neglected any 

variation of / caused by accretion. If M = M{t), 

we have z>(t) = (27r/)~i £o Afo (M(^)/Mo)^~"/^ where 

io = (GMi?mo)^^^, and i?mO and Mq are R^ and M at 
t = Tq, respectively. For the a — case we assumed 
^0 = ^max- In this case we therefore assume that the 
system is accreting the maximum specific angular mo- 
mentum possible, giving an upper limit on the spin-up 
torque and therefore a lower limit to the accretion rate 

Mo. 

Since we assumed M{t) oc L{t), to determine the tem- 
poral dependence of M{t) we studied the energy spec- 
tra of the source for each continuous interval of data 
combining PCA and HEXTE data. All the spectra 
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are well fitted with a model consisting of a power law 
with an exponential cutofF plus thermal emission from 
a Keplerian accretion disk modified by photoelectric ab- 
sorption and a Gaussian iron line. In order to derive 
L{t) for each spectrum we made the simple assumption 
L{t) oc -F'(3-i5o)(i), which is the unabsorbed flux in the 
RXTE PCA plus HEXTE energy band (3 - 150) keV. A 
good fit of -F(3-i5o)(0 vs t between December 7 and 14 
(Atobs ~ 7 days) is given by the expression -F(3_i5o)(t) = 
-P(3-i50) [1- (t-To)/tij] with tB = 8.4±0.1 days, where 
F(3_i5o) is the unabsorbed flux at t = Tq. Therefore 

we have v{t) = i>o [I - {t - ro)/^^]^""/^, where the spin 
frequency derivative at t = Tq is vq = (27r/)~^ £o Mq. 

With this expression for ^{t) eq. (2) can be readily 
integrated. Since e = (t — To)/tB < 1 for t < 7 days, 
we took a series expansion of the integral and obtained 

SMt) = l/2z>o(t-ro)2[l- (2-a)(f-ro)/(6tB)+e], 
with an error ^ < a(l — Q:/2)/24e^. 

We have therefore fitted these phases with the func- 
tion: 



through the formula: 
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Using the best fit value for Ai'o we computed the im- 
proved spin frequency estimate and repeated the same 
procedure described at the beginning of this paragraph, 
folding at the new estimate of the spin period. The new 
phases were fitted with eq. (6). In this case, At'o was fully 
compatible with zero. These phases are plotted versus 
time in Figure 1 (upper panel) together with the residu- 
als in units of a with respect to eq. (6) (lower panel). The 
best fit estimates of uq and !>o are reported in Table 1 for 
three values of a, namely a = which correspond to a 
location of i?ni independent of the accretion rate (cfr. the 
model of Rappaport Fregeau, and Spruit 2004 in which 
= ^co for any M), the standard case a = 2/7 which 
corresponds to i?m proportional to the Alfven radius, and 
a = 2 which has been given for comparison purposes and 
corresponds to a parabolic trend, expected in the case of 
constant M. Of course, the value of obtained in this 
latter case is in agreement with the value obtained by 
F05 of !> = 8.4(6) X 10~^3 Hz s~^ A comparison of the 
X^/dof for each of the adopted values of a (also reported 
in Table 1) shows that the statistics is not good enough 
to distinguish between these three possibilities. 

Finally, to evaluate the systematic errors on the spin 
frequency and its derivative, we adopted the positional 
uncertainty of 0.06" radius reported by Rupen et al. 
(2004) in a series expansion of eq. (5), finding aj/syst ~ 
2.2x 10~^ Hz and CT,>syst ~ 4.4 xlO"-^^ Hz/s, respectively. 
Even adopting a positional error circle of 0.2" (which is 
the distance between the optical and radio position, Fox 
& Kulkarni 2004), the resulting systematic uncertainties 
are (Jvsyst ~ 7.3 x 10~^ Hz and ai-syst ~ 1-5 x 10~^^ 
Hz/s. Note that the systematic error on the spin fre- 
quency is comparable with the error derived from the fit 
of the phase delays and reported in Table 1, the system- 
atic error on the spin frequency derivative is at least one 
order of magnitude below the error derived from the fit. 

4. DISCUSSION 

From the best-fit value of the spin frequency derivative 
z>o we can compute the mass accretion rate at t = Tq 
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5.9 X i>_i3 /45 m-''/^{Rco/Rrnof'^, (7) 



where M_io is Afo in units of 10^"'^" Mq yr~^, v-iz is vq 
in units of 10^^'^ s^^, and 145 is / in units of 10^^ g cm^. 
In the following we will adopt the FPS equation of state 
for the neutron star matter for m = lA and the spin fre- 
quency of IGR J00291-^5934 which gives 145 = 1.29 and 
R = 1.14 X 10^ cm (see e.g. Cook, Shapiro & Teukolsky 
1994). This gives a lower limit in the mass accretion rate 
of Af_io ^ 70 ± 10 (case a = 0). In order to compare 
the experimental estimate of Mq with the observed X-ray 
luminosity, we have to derive the bolometric luminosity 
L{t) from the observed flux J^(3_i5o)(i). To this end we 
consider the spectral shape at t = To in more detail. 

Since the value of the hydrogen column Nn is poorly 
constrained in the RXTE energy band which starts at 2.5 
keV, we fixed it to TVr = 2.80 x lO^^ cm'^, which is the 
value obtained by Nowak et al. (2004) analysing Chandra 
data. Our spectral results are practically independent 
P^f the precise value of the iVn below the total Galactic 
drogen column in the direction of IGR J00291-I-5934. 
The power law is the dominant spectral component. 
In particular we found a power law spectral index a = 
-0.59t^;^3^ and an e-folding energy Efoid ~ 178t^? keV 
(with values ranging from 60 to over 260 keV through- 
out our data). The ratio of the unabsorbed fluxes in 
the bands 0.001 - 1000 keV and 3 - 150 keV is 1.6. 
This ratio is almost independent on the e-folding en- 
ergy, increasing up to 8% when the e-folding energy 
increases from 60 to 260 kcV. Therefore we assume 
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The power-law component 



presumably originates in regions of small optical depth 
just above each polar cap (see e.g. Poutanen & Gierlinski 
2003; Gierlinski & Poutanen 2005), thus we neglect, to 
first order, any effect of the inclination of the emitting 
region with respect to the observer. On the other hand, 
we observe a single-peaked pulse profile, which means 
that we only see the emission from one of these regions 
(e.g. Kulkarni & Romanova 2005). 

If this is the case, we have to take into account the pos- 
sibility that we are underestimating the total flux in the 
power-law component because one of the two polar caps 
is never visible (this can happen if the sum of the angles 
between the magnetic axis and spin axis and between the 
line of sight and the spin axis is less than 7r/2). If we only 
see polar cap A, we indicate with Fa the flux emitted by 
A reaching the observer. Because the emission from A is 
isotropic, to obtain the luminosity we have to integrate 
over 47rd^, so we obtain La = FA4:TT(f. However, under 
our assumptions, this is an underestimate of the total (A 
+ B) luminosity of the system consisting of the two polar 
caps, since part of the flux (Fb) that should be emitted 
in the direction of the observer is never visible. 

Actually almost half of the emission from cap A fac- 
ing the neutron star surface is intercepted by the neu- 
tron star itself and re-emitted towards the observer, and 
this re-emission should indeed be consider. However, it 
is reasonable to assume that this intercepted emission 
is reprocessed by the neutron star and re-emitted as a 
blackbody-like spectrum at a relatively low temperature. 
This reprocessed emission has therefore a very different 
spectral shape from the original power-law component: 
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simple estimates demonstrate that the temperature asso- 
ciated with this blackbody-hkc emission is below 1 keV, 
and thus not related to the power-law component. More- 
over, most of the emission from this component is outside 
the energy range of RXTE/PCA (and thus poorly con- 
strained by the RXTE observation). 

Since the emission from an optically thin region is pro- 
portional to the volume of the emitting region visible, 
in the hypothesis that we totally miss the flux from one 
of the two polar caps, we should multiply by a factor 2 
the unabsorbed flux of the power law in order to take 
into account the emission of the optically thin region 
above the unseen polar cap. We therefore parametrize 
the luminosity of the power-law component with a fac- 
tor T] (which can assume values between 1 and 2), and 
we can write the total luminosity of the two polar caps 
as: LpL ^ r/FpL (o,oo) x ^^d^ = 0.75^°;?° x lO^^ rj d^^^^ 
erg/s, where rfskpc is the source distance in units of 5 
kpc. The uncertainty on the luminosity has been evalu- 
ated conservatively by propagating the uncertainties on 
the spectral parameters treated as they were independent 
of each other. 

The second component is the thermal emission inter- 
preted as emission from a Shakura-Sunyaev accretion 
disk, that is fitted with the diskbb model. We foimd 
a temperature of the inner rim of the accretion disc of 
Tin = 0.68^0 25 keV. On the other hand, because of the 
poor coverage of RXTE/PCA at soft X-rays, the disk 
blackbody normalization, K = (i?inkm/£'iokpc)^ cos^ = 
25^20'^ (where Diokpc is the distance in units of 10 kpc), 
is basically unconstrained. We therefore use the inner 
disk temperature and the Virial theorem to infer the 
bolometric luminosity of the disk, as follows. We in- 
terpret the inner temperature of the accretion disk as 
derived by the diskbb model as the maximum temper- 
ature in the disk. To prove this we have also fitted the 
soft X-ray spectrum with the diskpn model (instead of 
the diskbb model): diskpn takes into account correc- 
tions for temperature distribution near the compact ob- 
ject using the Paczynski-Wiita pseudo- Newtonian poten- 
tial (sec Gicrlinski et al. 1999). Also this model gives a 
maximum temperature of the disk of 0.70 keV (fully com- 
patible with the value obtained with the diskbb model). 
Standard disc theory (see e.g. Frank, King, & Reine 
2002) predicts that the temperature of the disk attains 
a maximum value at a radius of 49/36 i?, , where i?, is 
the radius at which the disk is truncated, correspond- 
ing to RmO in our case. Using this in eq. 5.43 of Frank, 
King, & Reine (2002) we obtain the inner disc radius 

through the relation: R^o = 1.77 m^/^ M]_{lT~^y^^ km. 
Combining this with eq. 7, and adopting Tkev = 0.70, 
we solve for the mass accretion rate and the inner disc 
radius. We foimd a mass accretion rate M-w = 85 ± 19 
and an inner disc radius i?mO — l-46tQ!4g ^ cm, 
that this is exactly in the very narrow range between 
the neutron star radius (~ 10^ cm) and the corota- 
tion radius (~ 2.4 x 10^ cm); the agreement with the 
expectation is compelling. The Virial theorem allows 
to calculate the fraction of the total luminosity that 
is emitted by the disc: 0.5i?/i?niO = 0.39. Therefore 
Lbbo = 0.39/(1 - 0.39) x Lplo = 4.8 x lO^S/yd^^p^ 
erg/s. 

The total bolometric luminosity is therefore Lq = 



1.23;o;i^ x lO'^^ r/rf^j^p^ erg/s. If we compare this lu- 
minosity with the mass accretion rate inferred from the 
timing analysis (assuming an efficiency C = 1) we obtain 
a distance to the source in the range (10.5-15)/?7i/2 kpc. 
Note that 10 kpc is close to the edge of our Galaxy in 
the direction of ICR J00291-h5934. If we push the factor 
rj to its maximum value of 2, we obtain a more reason- 
able range of distances to the source of 7.4 — 10.7 kpc, 
consistent with the lower limit of 5 kpc discussed in F05. 
Note also that the effect of including magnetic torques 
due to threading (see assumption iv in § 3) can only push 
the source further away. This is evident from the a = 
case (which gives a spin derivative similar to the a = 2/7 
case discussed here), in which we assume that the sys- 
tem is accreting the specific angular momentum at the 
corotation radius, that is the maximum possible; since 
any torque beyond the corotation could only spin down 
the system, this would increase the required M to justify 
the measured spin-up. 

Finally, using the value of the inner disc radius de- 
rived above and the normalization of the disc blackbody 
model, K = (-Rinkm/-Diokpc)^ cos = 25^20°' i^^^^r 
the inclination of the system with respect to the line of 
sight, which is i > 40° for a distance of 9.5 kpc. 

5. CONCLUSIONS 

We have analysed RXTE data of the fastest known 
accreting millisecond pulsar, ICR J00291+5934, during 
the period 7 — 14 December, 2004. We report a revised 
estimate of the spin period and discuss the spin period 
derivative. The source shows a strong spin-up, which 
indicates a mass accretion rate of about 8.5 x 10~^ Mq 
yr~^. We have checked that this mass accretion rate 
is compatible with the X-ray spectrum of the source. 
In particular we have shown that, with this high accre- 
tion rate, the measured temperature of the disc black- 
body emission implies an inner disc radius in excellent 
agreement with the inferred magnctosphcric radius of the 
source (constrained to be in the narrow range (1 — 2) x 10^ 
cm) and that the source is probably seen at high incli- 
nation. However, the mass accretion rate inferred by the 
observed spin-up (and calculated using standard values 
for the neutron star moment of inertia) , using the simple 
relation Lboi = C,GMM / R (with an efficiency ^ = 1), 
would correspond to a quite high bolometric luminosity 
of the source of Lboi ~ 10'^* ergs/s, much higher than 
the observed source luminosity assuming a distance of 5 
kpc. We propose that the simplest explanation of this 
discrepancy is that part of the accretion luminosity is 
not visible. Indeed, if we only sec the emission of one 
of the two polar caps, we could miss up to half of the 
flux in the power-law component. Under this hypothesis 
(described in detail in the previous section) , we have ex- 
trapolated the flux of the source derived from the X-ray 
spectrum, which corresponds to a bolometric luminosity 
ofLboi — 2.5xl0''''(i^]^p^ erg/s. Comparing this extrapo- 
lated bolometric hmiinosity with the mass accretion rate 
of the som'ce as derived from the timing, we find an agree- 
ment between these two quantities if we place the source 
at a distance between 7 and 10 kpc. 

Other possibilities to explain this discrepancy can be 
that part of the accretion luminosity is not observed be- 
cause emitted in other energy bands or because the effi- 
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Table 1. Orbital and spin parameters of IGR J00291+5934. 



G05 



F05 



This work 



oi sin i (It-ms) 

(MJD) 



Projected semimajor axis 
Orbital period, -Porb (s) 
Epoch of ascending node passage, T* 
Eccentricity, e 
Spin frequency, uq (Hz) 
Spin frequency derivative 
Spin frequency derivative 
Spin frequency derivative 
Epoch of the spin period. 



i>0 (Hz/s) (i/ = constant) 

I/O (Hz/s) (a = 0)'' 
i>0 (Hz/s) (q: = 2/7)*^ 
To (MJD) 



64.993(2) 

8844.092(6) 

53345.1619258(4) 

< 2 X 10"* (3 a) 
598.89213064(1) 

< 8 X 10-13 (3 cr) 



598.89213060(1) 
8.4(6) X 10"" 



53346.0 



598.89213053(2) 
0.85(0.11) X 10-12 
1.17(0.16) X 10-12 
1.11(0.16) X 10-12 
53.346.184635 



(x2/dof=106/77 ) 
(x2/dof=113/77 ) 
(x2/dof=lll/77 ) 



Note. — Errors arc given at la confidence level. The errors quoted for the spin frequency and spin frequency derivative are derived from the 
phase delays fitting, and do not include systematic errors induced by the source position uncertainty. Adopting a positional uncertainty of 0.06" 
(Rupen et al, 2004), these are: <t„ syst ~ 2.2 X 10^** Hz and a,jsyst ~ 4.4 X lO^^'^ Hz/s, respectively. 

^G05 reported a value of 53345.1875164(4) MJD for the epoch of superior conjunction, i.e. when tiie NS is behind the companion: as, 
in this work, we considered the epoch of ascending node passage as a reference time, the G05 reference time reported here has been 
decremented by Porb/4, 

''Averaging over the 7 days of our observation, we get < i>o >— 0.68 x IQ-i^ Hz/s. 



ciency ( of the conversion of the gravitational potential 
energy of the accreting matter into X-ray luminosity is 
less than 1, or because of occultation effects (which may 
be favoured if indeed the source is highly inclined). In 
these cases wc should conclude that the observed X-ray 
luminosity is not a good tracer of the total mass accretion 



rate, M, onto the neutron star. 
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Fig. 1. — Pulse phases computed folding at the spin period reported in Table 1 and plotted versus time together with the best fit curves 
(upper panel) and residuals in units of cr with respect to the model with a = 2/7 (lower panel). 



